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We study the properties of coupled linear and nonlinear resonances. The fundamental phenomena 
and the level crossing scenarios are introduced for a nonlinear two-level system with one decaying 
state, describing the dynamics of a Bose-Einstein condensate in a mean-field approximation (Gross- 
Pitaevskii or nonlinear Schrodinger equation). An important application of the discussed concepts 
is the dynamics of a condensate in tilted optical lattices. In particular the properties of resonance 
eigenstates in double-periodic lattices are discussed, in the linear case as well as within mean-field 
theory. The decay is strongly altered, if an additional period-doubled lattice is introduced. Our 
analytic study is supported by numerical computations of nonlinear resonance states, and future 
applications of our findings for experiments with ultracold atoms are discussed. 
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I. INTRODUCTION 

In the last decade, the advance of atom and quantum 
optics has made it possible to realize and to study the 
evolution of the center-of-mass motion on scales rang- 
ing from the microscopic (single particle) to the macro- 
scopic (many-particle) realm [ij, |2j . In a typical exper- 
iment with ultracold atoms, interactions can either be 
made negligibly small or reduced to a mean-field effect 
on the evolution of the macroscopic order parameter of 
a Bose-Einstein condensate (see, e.g., 0, Q and refer- 
ences therein). The latter approach results in an effec- 
tive nonlinear Schrodinger equation, the following Gross- 
Pitaevskii equation: 

(i) 

that describes the dynamics of the macroscopic wave 
function (or of the order parameter) of a Bose-Einstein 
condensate (BEC) for zero temperature Q. This mean- 
field description has proved to be extremely successful 
and reliable for most of recent experiments. The non- 
linearity of the equation leads to a variety of surprising 
phenomena, which are present even in a simple nonlinear 
two-level system. Self-trapping of a BEC in a double-well 
trap was observed experimentally only recently |4] . The 
self-trapping transition manifests itself in the appearance 
of novel nonlinear eigenstates Q. The appearance and 
disappearance of nonlinear eigenstates may also lead to a 
breakdown of adiabaticity and nonlinear Zener tunneling 

sag. 
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In the present paper, we investigate nonlinear quantum 
dynamics in decaying systems. Up to now, only relatively 
few papers have studied nonlinear and non-hermitian 
quantum dynamics, discussing self-stabilizing, shi fting 
and broadening of nonlinear resonances |9L flOt fill fl3 . f]~^ | - 
Here we focus on the coupling of nonlinear resonances in 
nonlinear, non-hermitian level crossing scenarios. Our 
first object of investigation, the nonlinear two-level sys- 
tem with one decaying level, offers analytic access to this 
subject. The eigenvalues and eigenvalues of its linear 
counterpart show some interesting features, such as ex- 
ceptional crossing scenarios |l4| . 

A very natural experimental setup leading to nonlinear 
dynamics and decay is the dynamics of a Bose-Einstein 
condensate in a tilted or accelerated optical lattice, cor- 
responding to the Wannier-Stark scenario of solid-state 
physics [lq . 

The decay dynamics in a nonlinear Wannier-Stark sys- 
tem was recently discussed in 0, It was shown 
that a nonlinear mean-field interaction can destroy res- 
onant tunneling. In this paper, we extend these studies 
to a double-periodic optical lattice. The decay dynam- 
ics in this system shows some interesting features even 
in the linear case, such as a splitting of resonant tunnel- 
ing peaks. The different types of non-hermitian crossing 
scenarios can be observed in dependence of the system 
parameters. 

The paper is organized as follows: first of all we re- 
view some important results about the crossing scenar- 
ios in the non-hermitian two-level system in the linear 
(section II) and the nonlinear (section III) case. The 
double-periodic Wannier-Stark system is introduced and 
analyzed in section IIVI Nonlinear Wannier-Stark reso- 
nances for a doubly periodic lattice are presented in sec- 
tion A discussion of novel experimental applications 
of our findings follows in section VI. 
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FIG. 1: Real (left) and imaginary (right) part of the eigen- 
values as a function of e for 7 = 1. A type I crossing is 
found for v = 1.01 > 7 (upper figures), a type II crossing is 
found for v — 0.99 < 7 (lower figures). 
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FIG. 2: Real (left) and imaginary (right) part of the eigen- 
values of the four-level Hamiltonian Q as a function of e for 
j = l,v = 0.98, w = 0.04 and 6 = 0.2. 



II. CROSSING SCENARIOS OF RESONANCES 
IN LINEAR QUANTUM MECHANICS 

We prepare for the full discussion of nonlinear reso- 
nance states as solutions of Eq. (|TJ by reviewing some 
essential properties of the simpler linear case. First of all, 
we want to illustrate the different types of possible curve 
crossing scenarios for non-hermitian systems. To start 
with, we briefly review a simple and instructive model 
system, a two-level Hamiltonian with one decaying level 

H3E3 

with e, v, 7 £ K and 7 > 0. In this approach it is assumed 
that one of the bare states decays with rate 7, while the 
decay is negligible for the other one. The mean energy 
of the two bare states is set to zero, the energy difference 
is given by 2e. The two states are coupled with strength 
v. A different, non-hermitian two-level Hamiltonian was 
previously discussed by Berry [l7j . 

The eigenvalues of the non-hermitian Hamiltonian 
are given by 

£± = -i 7 ± y/(e-i-y) 2 + v 2 =E±- iT ± /2. (3) 

Both real and imaginary part of the eigenvalues are differ- 
ent for e =/= 0. (Anti) crossings of the real and imaginary 
part are found only in the critical plane e = 0. The ex- 
ceptional line v = ±7 separates the critical plane into 
different regions: 

(a) For |u| > 7 the imaginary parts of the eigenvalues 
coincide, T + = T_ = 27, while the real parts differ. 
This case is denoted as type I crossing. 

(b) For |u| < 7 the real parts of the eigenvalues co- 
incide, E+ = E— = 0, while the imaginary parts 
differ. This case is denoted as type II crossing. 

(c) The eigenvalues are fully degenerate, £ + — £-, 
along the critical lines v — ±7. 



The two different crossing types are illustrated in Fig. ^ 
For a type I crossing, i.e. \v\ > 7, the imaginary parts of 
the eigenvalues cross while the real parts anti-cross. For 
a type II crossing, i.e. |u| < 7, it is the other way round. 
Physically this crossing describes a resonantly enhanced 
tunneling (RET) effect: the decay rate of the lower state 
increases significantly if this state is energetically close 
or equal to the decaying upper level. 

In view of the discussion of Wannier-Stark resonances 
in period-doubled lattices in the following sections we 
want to introduce another model system. We assume 
that the bare states split up into two states, where the 
energies of the stable bare states differ slightly by 25. 
Each stable state mainly couples to one of the decaying 
states, while all other couplings are assumed to be weak. 
We consider the Hamiltonian 

n _(H 2 +A W \ , , 

Hi ~{ W H 2 -Aj ^ 

with the two-level Hamiltonian H2 defined in (J2J) and 

*-'(S:)- "-(!!)■ (5) 

Figure [5] shows the eigenvalues of this Hamiltonian in de- 
pendence of the on-site energies e. One observes that the 
resonance peak splits up into two peaks. Two possibil- 
ities for resonant tunneling, i.e. two type II crossings, 
are found instead of just one. This crossing scenario is 
robust against small variations of the coupling w as long 
as w <C v is fulfilled. A non- vanishing coupling w causes 
as slight asymmetry of the two crossings. For v > 7 one 
has two type I crossings instead, i.e., the imaginary parts 
of the eigenvalues cross while the real parts anti-cross. 

The change of a system parameter, e.g. the strength of 
the Stark field F in the Wannier-Stark system discussed 
in section llVl will typically affect the bare state energies 
e as well as the decay rate 7 and the coupling strengths. 
Therefore we consider a variety of the four-level Hamil- 
tonian PJl. the parameters of which are functions of the 
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e = —F/2 + 0.1 FIG. 4: Real (left) and imaginary (right) part of the eigen- 

values as a function of e for 7 = 1.00 and c = 0.3 and 
1 — F c v _ ^ Qcj ^ ^ ( U pp er figures), v = 1.00 = 7 (middle figures), 

v = 0.05 F c~ 1 / 2F v = 0.95 < 7 (lower figures). The linear levels (c = 0) are 
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— 01 F e~ 1 ^ 2F (6) plotted as dotted lines for comparison. 



The exponential scaling of the decay rate 7 is well known 
from standard Landau-Zener theory (see, e.g., 01)- The 
dependence of the bare state energies and the coupling 
coefficients on F were analyzed in detail for a two-ladder 
system in [l^. It was shown that the Wannier-Stark 
spectrum is accurately described assuming that the pa- 
rameters scale as in Eq. ©. The actual values of co- 
efficients in © are chosen in an ad hoc manner for il- 
lustration only. The resulting decay rates are illustrated 
in Fig. |2| For S = 0.02, we find two type II crossings. 
With increasing S, the crossing on the right changes its 
form and becomes a type I crossing. Such a crossing sce- 
nario is naturally realized for Wannier-Stark resonances 
in double-periodic lattices, as will be shown in Fig. |5] be- 
low. 



III. NONLINEAR NON-HERMITIAN 
CROSSING SCENARIOS 

The linear non-hermitian two-level system described 
in section ITT1 neglects any particle-interaction. Including 
this interaction in a meanfield-description according to 
the Gross-Pitaevskii-equation Q yields a nonlinear non- 
hermitian two- level system [2(| , described by the Hamil- 
tonian 

( e + 2c|^ 1 | 2 -2i 7 v \ 

where the nonlinearity parameter c is proportional to 
the parameter g in the Gross-Pitaevskii-equation l[T]l. 
The nonlinear eigenstates are then defined as the self- 



consistent solutions of the time-independent Gross- 
Pitaevskii-equation 

«-(£)-"(£)• < 8 » 

The nonlinear eigenstates crucially depend on the nor- 
malization of the state vector, which is fixed as "01 1 2 + 
\ipi\ = 1 throughout this section. For convenience we 
symmetrize the nonlinear Hamiltonian (JJJ by substract- 
ing a constant energy term c(|-0i| 2 + IV^I 2 )- The Gross- 
Pitaevskii-equation |JSJ then reads 

( e+ T 2i7 _:„)(£)-"(£) w 

with k = \ipi\ 2 — l^l 2 - The self-consistent solutions of 
this nonlinear equation define the nonlinear eigenstates 
and eigenvalues. 

Note that the nonlinear eigenstates are not connected 
to stationary solutions of the time dependent system, if 
the chemical potential turns out to be complex, since the 
dynamics depends crucially on the normalization of the 
state- vector, which is not constant for a complex valued 
chemical potential. After some algebraic manipulation 
one can show that the nonlinear eigenstates, i.e., the so- 
lutions of equation @ are given by the real roots of the 
equation 

(c 2 + 7 2 )k 4 + 2cen 3 + (v 2 + e 2 - 7 2 - c 2 )n 2 - 2cen - e 2 = 0. 

(10) 

Depending on the parameters, there are two or four real 
roots and each of them is connected to a complex eigen- 
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FIG. 5: Real (left) and imaginary (right) part of the eigen- 
values JOJ as a function of e for 7 = 1.00 and c = 0.9 and 
v = 1.01 > 7 (upper figures), v = 1.00 = 7 (middle figures), 
v = 0.99 < 7 (lower figures). 



value by 



H = c + c/k- + k) = M - ir/2. 



(11) 



For e — > ±00, the linear term dominates and one has only 
two eigenvalues. For e = 0, K = is a double degener- 
ate solution of equation I|1U|) . For 7 < v these states are 
connected to the common linear (anti-)symmetric eigen- 
states, while this is not the case for 7 > v. In the follow- 
ing we consider the crossing scenario of the eigenvalues 
in dependence on e for different fixed values of the other 
parameters. 

The nonlinear eigenstates of a two- level system are well 
known for the hermitian case 7 = [1 H El. Novel 
eigenstates emerge with broken symmetry if the nonlin- 
earity exceeds a critical value, |c| > c cr = u, which is 
given by the coupling strength v which corresponds to 
half of the gap between the linear levels at e = 0. The 
levels show looped structures around e = with a width 

|e| < (c 2 /3„ ,,2/3)3/2^ 

Let us first discuss the effect of a weak nonlinearity: 
Figure 01 shows the eigenvalues for a relatively weak non- 
linearity in comparison with the linear case c = 0. The 
size of the gaps is not altered by the nonlinearity, which 
leads to the important fact that the nonlinearity does 
not influence the crossing type. Nevertheless it changes 
the shape of the levels. For a type I crossing, the real 
part of the upper level is sharpened while the one of the 
lower level is flattened, which is well-known for 7 = 
;6| . For a type II crossing the effect is basically the same, 
but is accompanied by a shift of the crossing point from 
jit = to fi = c and the lower level is stretched to this 
point. The imaginary parts bend slightly to the left. At 
the exceptional point additional eigenvalues emerge in a 



narrow interval around e = 0. 

In general, the presence of moderate decay facilitates 
the formation of novel eigenstates. Figure [S] shows the 
nonlinear levels for c = 0.9, which is slightly below the 
critical value c cr = v for 7 = 0. For a type I crossing, 
v > 7, the real parts show a familiar loop. The imagi- 
nary parts cross as usual, showing an additional S-shaped 
structure. In fact, the critical nonlinearity for the emer- 
gence of looped levels is decreased to 



\Jv 2 - 7 2 , 



(12) 



which can be seen by analyzing the behavior of the poly- 
nomial (|10|) for e = 0. At the exceptional point v = 7, the 
critical nonlinearity tends to zero and there are additional 
eigenstates even in the case of arbitrary weak nonlinear- 
ity. For a type II crossing, v < 7, the imaginary parts 
anti-cross and the real parts cross in a manner which can 
be understood as a destruction of the loop at its lower 
edge. The crossing appears at the former crossing point 
of the loop at e = and /i = c 0. In the non-decaying 
case, 7 = 0, novel eigenstates firstly emerge at the point 
of the avoided crossing at e = 0. This remains true for 
a type I crossing, v > 7. For a type II crossing, v < 7, 
however, novel eigenvalues emerge, again in an S-shaped 
structure around some non-zero value of e. 

Concluding this section, a weak nonlinearity does not 
alter the crossing type, however it deforms the levels in 
a characteristic manner. For a type II crossing the real 
parts cross at /1 = c ^ 0. At the exceptional point novel 
eigenstates emerge, even for small nonlinearities. The 
presence of decay facilitates the formation of novel eigen- 
states for stronger nonlinearities. For type I crossings, 
loops appear if |c| > c cr = y/v 2 — j 2 . For a type II 
crossing the additional eigenstates emerge around some 
nonzero value of e forming a double S-structure. If the 
sign of the nonlinearity is changed, the levels interchange 
their behavior, i.e., the real parts are mirrored at the 
e-axis, the imaginary parts at the T-axis. 



IV. WANNIER-STARK RESONANCES IN 
DOUBLE-PERIODIC LATTICES 

A. Fundamentals of the linear Wannier-Stark 

system 

A prime example for resonances and resonant tunnel- 
ing is the (linear) Wannier-Stark problem described by 
the Hamiltonian 



//u s = ~\ih +v{x)+Fx 



(13) 



with a periodic potential V(x + d) = V(x). We use 
rescaled units in which % = M = 1. The Wannier-Stark 
problem was already discussed in the early days of quan- 
tum mechanics in the context of electrons in solids under 
the influence of an external electric field j^]. Coher- 
ent dynamics of electrons in semiconductor superlattices 



5 




FIG. 6: Eigenenergies (upper panel) and decay rates (lower 
panel) of the two most stable Wannier-Stark ladders for the 
potential V(x) — cosx. 

were observed not until the 1990's. Experiments showed 
Bloch oscillations for 'weak' electric fields and decay for 
stronger fields The Wannier-Stark system is fur- 

thermore realized for the propagation of light pulses in 
thermo-optically biased coupled waveguides. Bloch oscil- 
lations as well as decay could thus be observed directly 
in real space |25j | . On the other hand, recent experiments 
with cold atoms and Bose-Einstein condensates in optical 
lattices offer some considerable advantages [21 1261 127^ l28j . 
Scattering by lattice defects or impurities is absent and 
experimental parameters can be tuned in a wide range. 
The periodic potential is generated by a standing laser 
beam and thus simply cosine shaped. 

Let us briefly review some fundamentals of Wannier- 
Stark resonances, which are defined by the eigenvalue 
equation 

H W S^ a A X ) =£a,n*a,n(x). (14) 

Here, a is the ladder index and n £ Z labels the lat- 
tice sites. The Wannier-Stark Hamiltonian 1)13(1 is non- 
hermitian due to the boundary condition: A wave packet 
will eventually decay towards x — > — oo. In fact, it has 
been shown that the spectrum of the Hamiltonian is con- 
tinuous with embedded resonances [2j|. Thus the reso- 
nance eigenenergies are complex, £ a<n — E a>n — \Y a /2, 
where the imaginary part Y gives the decay rate. The 
Wannier-Stark Hamiltonian has one important symme- 
try, it is invariant under a simultaneous spatial transla- 
tion over a lattice period d and an energy shift dF. This 
symmetry is expressed by the commutation relation 

[Hws,T m ] = -mdFT m , (15) 

where T m is the translation operator over m lattice peri- 
ods. Now it is easy to see that the Wannier-Stark reso- 
nances from one ladder a are related by a simple trans- 



lation, 

H W sT m ^ a ,n{x) 

= T m H ws ^ a .n{x) + [Hws,T m ]^ a ,„(x) 

= {£ a}n -mdF)T m V atn (x). (16) 

Thus the discrete spectrum is arranged in the form of the 
so-called Wannier-Stark ladders, 

£ a , n = E a , + ndF - iY a /2. (17) 

The different ladders are labelled by a = 0,1, 2,.... An 
efficient method to calculate the resonance eigenstates 
was introduced in |3p| . a recent review can be found in 

m 

The decay rate Y a is the same for all resonances in 
one ladder. In general, the decay rate scales as Y ~ 
F exp(— irAE 2 /F), where AE is the energy gap between 
the Bloch bands of the periodic potential. This result can 
be deduced from Landau-Zener theory 0, [l^. How- 
ever, one observes peaks of the decay rate on top due 
to resonant tunneling. For example, the decay rate of 
the two most stable resonances for the periodic poten- 
tial V(x) = cosx is plotted as a function of the inverse 
field strength 1 jF in figure HO RET takes place when 
a state of a lower ladder with energy £ Q) „ gets in reso- 
nance with a state of a higher ladder at a different site, 
i.e. E a ^ n = E a i. n i. The decay rate of the lower ladder is 
significantly increased as it couples resonantly to a higher 
ladder with a higher decay rate. For example, the pro- 
nounced peak in the ground ladder decay rate (i.e. a = 0) 
at F « 1/7 corresponds to the resonance a = <-» a' = 1 
and n' = n — 1 . 

B. Double-periodic lattices 

Now we turn to the main subject of the present paper. 
We consider a double-periodic potential V(x) consisting 
of a major optical lattice of period d plus an additional 
shallow lattice of double period, 

V{x) = V (sm 2 (nx/d) + <5 sin 2 (7ra;/2d + </>/2)) . (18) 

Rescaling the spatial coordinate as x' — 2irx/d and ne- 
glecting a constant potential offset, we can rewrite the 
periodic potential as 

Vn 

V{x) = -y (cos(x) + 5cos(a;/2 + </>)) . (19) 

The relative phase of the two lattices is denoted by <f>. 

Due to the additional lattice each Bloch band splits up 
into two minibands , and each Wannier-Stark ladder 
splits up into two miniladders, as proved in the follow- 
ing. The symmetries of the Hamiltonian are given by the 
commutation relations 

[H, T 2m ] = -2mdFT 2m 

[H, T 2m+1 G] = -(2m + l)dFT 2m+1 , (20) 
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FIG. 7: (Color online) Resonant tunneling in a double- 
periodic Wannier-Stark system. Decay rates of the four most 
stable resonance for Vo — 2, <j> — and 5 = (dashed black 
line), 5 = 0.05 (solid red lines) and S = 0.1 (solid blue lines). 

where T m is the translation operator over m lattice peri- 
ods and the operator G inverts the sign of 5 in all follow- 
ing terms. Then it is easy to see that the Wannier-Stark 
states of one ladder are related by a translation over an 
even number of lattice periods, or by a translation over 
an odd number of lattice periods plus an inversion of the 
sign of 6, 

= (£ a , n (6) - 2mdF) T 2m V a , n {x), (21) 
HT 2m+1 G^ a , n (x) 

= (£ a ,n(S) - (2m + l)dF) T 2m+1 G* Q ,„(x). 

Furthermore, it can be shown that the energy offset £ Qi o 
is antisymmetric in 5, 

£ a ,o{- S ) = -£*A S )- ( 22 ) 

Thus the Wannier-Stark ladders split up into two mini- 
ladders, each with an energy offset 2£ a (6): 

£a,2n = £a{5) + 2ndF 

£ a an+i = -£ a (S) + (2n+l)dF. (23) 

A similar proof is given in [32j within the tight-binding 
approximation. 

The decay of the Wannier-Stark resonances is seriously 
influenced by the additional period-doubled potential. 
Figure [7| shows the decay rate T as a function of the 
inverse field strength 1/F for Vo — 2, <fi = and 5 = 0.05 
and 5 = 0.1, respectively. The decay rate of the single- 
periodic lattice S = is also plotted for comparison. As 
all Wannier-Stark ladders split up into two miniladders, 
so does the decay rate T. The general scaling of T with 
F remains the same for both miniladders, while the RET 
peaks are seriously altered. The peaks split up into two, 
where the height of the subpeaks increases significantly. 
This effect is mostly pronounced for the major resonance 
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FIG. 8: (Color online) Explanation of resonantly enhanced 
tunneling (RET). See text for details. 



at F w 1/7. The explanation of the splitting is straight- 
forward: In figure [SI the solid blue lines represent the 
energy levels of the two most stable Wannier-Stark lad- 
ders. Due to the period-doubled potential they are al- 
ternately shifted up or down with respect to the unper- 
turbed level (5 = 0, dash-dotted lines). The decay rate T 
is resonantly enhanced when two energy levels of differ- 
ent ladders are degenerate, as indicated by the solid red 
arrow. However, for 5 ^ 0, the two miniladders are in 
resonance with higher ladders for different values of the 
field strength F due to the alternating energy shift. In 
figure |HJ for example, one of the miniladders is off res- 
onance (indicated by the broken arrow) while the other 
one is in resonance (solid arrow). If the field strength is 
lowered, the off-resonant ladder will get into resonance 
at another value of F. 

With increasing amplitude 5 of the second double- 
periodic lattice, the splitting of the resonant tunneling 



5 = 0.075 < n -i 




FIG. 9: (Anti) Crossing of the real and the imaginary part 
of two eigenenergies of the Wannier-Stark resonances in a 
double-periodic Wannier-Stark system for Vo — 2, (f> = and 
5 — 0.075 (upper panels), S — 0.08 (lower panels). 






FIG. 10: (Color online) Upper panel: Decay rates of the two 
most stable resonances as a function of the relative phase <f> 
and the inverse field strength 1/F for Vo — 2 and 8 = 0.05. 
Lower panel: For fixed Stark field F = 1/9. 

peaks clearly becomes more pronounced as one can see 
in figure [3 Another (not so intuitive) effect is that the 
additional lattice can also alter the crossing type. For 
5 = 0.1 one observes a crossing of the decay rates at 
F w 1/9.5 due to resonant tunneling instead of an anti- 
crossing. This effect is further illustrated in figure^ The 
real and the imaginary part of the resonance eigenener- 
gies are plotted in the vicinity of one of the RET peaks 
(cf. figure EJ. For S = 0.075, one observes a familiar 
RET peak, i.e., a type II curve crossing. For 6 = 0.08, 
however, the crossing type is altered from type II to type 
I. The decay rate of one miniladder crosses the decay 
rate of one excited miniladder. Correspondingly, the real 
parts anti-cross. A diabolic point, where real and imag- 
inary part are degenerate, is found at S = 0.0772 and 
F = 1/8.937 for = 0. If we consider the relative phase 
cj) as another free parameter, the set of diabolic points is a 
one-dimensional subset of the three-dimensional parame- 
ter space (6, F, </>) ; the diabolic crossing has codimension 
one. 



C. Output control by the relative phase 

Up to now, we have shown that the RET peaks split 
up, whereby the splitting increases with the amplitude S 
of the additional lattice. Furthermore, the decay depends 
crucially on the relative phase <f> 01 the lattices. Figure 
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FIG. 11: Effects of fluctuations of the relative phase <f>. (a) 
Decay rate of the two most stable resonances for Vb = 2, 
6 — and <j> — 0. (b) Position Fi ( 2 of the RET peaks vs. 
the relative phase 4>. (c) Uncertainty AF of the positions of 
the RET peaks vs. fluctuation A(j> of the phase, (d) Relative 
fluctuation AF/To of the decay rate for F = 1/8 vs. Acj>. 



^| shows the decay rate of the two lowest miniladders in 
dependence of the relative phase <f> and the inverse field 
strength 1/F for V = 2 and S = 0.05. The splitting of 
the RET peaks is maximal for cf> = 0. It becomes zero for 
<f) w 7r/2, where the decay rates of the two miniladders 
degenerate again. Despite the fact that the additional 
lattice is much weaker than the single-periodic one (S = 
0.05), it seriously affects the decay properties. The decay 
rate of the two lowest miniladders vary over several orders 
of magnitude. This is further illustrated in the lower 
panel of figure fTOI where the decay rate is plotted for a 
fixed value of F = 1/9. The decay rate for the single- 
periodic case S — is also plotted for comparison. The 
decay rate as a function of the phase, r(0), is 7r-periodic. 
This can be seen as follows: Shifting the phase by an 
amount of — n and the position by 2-7T leaves the system 
invariant. A spatial translation by 2-7T exchanges the two 
miniladders. 

The sensitive dependence on the relative phase <f> might 
be very useful for the design and control of future exper- 
iments with double-periodic optical lattices. A complete 
stabilization of the relative phase of the two lattices will 
be hard to achieve and fluctuations will play an impor- 
tant role. Yet, Ritt et al. have realized a new technique 
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based on Fourier synthesis to control two (or even more) 
optical lattices with different spatial harmonics ^5^. In 
the following, we will discuss the effects on the decay 
rates of the Wannier-Stark resonance in some detail. The 
strongest effects are found in the vicinity of the RET 
peaks, on which we focus in our present discussion. 

Figure ITT1 fa) again presents the decay rate as a func- 
tion of the inverse field strength 1/F for Vo = 2 and a 
weak additional lattice with S = 0.05 and = 0. For 
the given parameters, we find RET peaks at F\ rs 1/6.9 
and i<2 w 1/8.4. Figure ITT1 (b) shows the position F% t 2 of 
the two the resonant tunneling peaks in dependence of 
the relative phase <fi. The positions vary in an interval of 
width AF w 1/40. 

In a real-life experiment, it is difficult to exactly control 
the relative phase of two independent standing waves. 
Therefore, we study also the influence of random phase 
fluctuations. If we assume that the phase <fi fluctuates 
in an interval of width A(f> around the desired value (f>o, 
4> & [4>q — A0/2 , </>o + A0/2] , the positions of the resonant 
tunneling peaks will also fluctuate in an interval of width 
AF. Fienre ITTl fc) shows how the width AF depends on 
the strength of the phase fluctuations A(j) for 4>o = and 
cf>o = 7r/2. The fluctuations are rather weak for <fio = 0, 
where the RET peaks have maximum distance. For a 
given value of the external field F, the decay rates T 
fluctuate in an interval of width Ar. FigurclTTlfd') shows 
the relative strength of the fluctuation, Ar /To vs. the 
phase fluctuations A<fi for the most stable resonance for 
an external field F — 1/8 and <fio — or <fio = tt/2, 
respectively. The relative uncertainty of the decay rate 
becomes greater than unity already for small fluctuations 
of the phase </>. Nevertheless, the noise induced shift and 
change in height of the RET peaks is small considering 
the absolute change of the decay rates around the RET 
peaks of about two orders of magnitude. 

However, it is also possible to exploit the sensitive de- 
pendence on the phase </>, if it can be accurately con- 
trolled. For instance, it could be possible to rapidly tune 
the output of a pulsed atom laser. An example will be 
discussed in detail in section IVll 



V. NONLINEAR WANNIER-STARK 
RESONANCES IN DOUBLE-PERIODIC 
OPTICAL LATTICES 

A method to obtain accurate, nonlinear Wannier-Stark 
resonances was proposed recently in 13], and we use a 
similar approach to numerically compute decay rates of 
the nonlinear version of the Wannier-Stark introduced in 
section llV Bl above. In contrast to the case studied in [l3| . 
the computations based on the Gross-Pitaevskii equation 
||TJ in the presence of a two-period optical lattice are 
more difficult, since the algorithm needs to discriminate 
between the two miniband solutions which are quite close 
in energy (cf. Figs. E]and0l. In particular, for very 
small Stark fields, it is hard to obtain convergence. In 
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FIG. 12: (Color online) Decay rates of nonlinear Wannier- 
Stark resonances in double periodic optical lattices with Vo = 
2, S = 0.05, 4> = and g = (dashed black line), g = -0.1 
(dashed green) g = —0.2 (dash-dotted red) and g = —0.3 
(solid blue), respectively. 



the following, we concentrate therefore on RET peaks 
at fields as large as possible, and one trick to shift the 
peaks to such values is to use attractive interactions, i.e. 
negative nonlinearities (g < 0). Figure IT2"1 presents a set 
of RET peaks for different values of the nonlinearity g. 

The nonlinearity induced a shift of both RET peaks 
(corresponding to the two minibands) and also a system- 
atic stabilization (i.e., smaller heights) can be observed, 
as predicted by similar results for the usual, one-band 
Wannier-Stark system [l^. This is analyzed in more de- 
tail in figure IT31 where the peak positions F Tes and the 
height of the peaks T los are plotted in dependence of g. 

The stabilization of the Wannier-stark states by an at- 
tractive nonlinearity is shown in the upper panel. An 
asymmetry of the two peaks is observed already in the 
linear case g — 0: The left peak is slightly higher than 
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FIG. 13: Shift of the RET peaks. Shown is the position F res 
of the RET peaks (lower panel) and the resonant decay rate 
r rcs (upper panel) in dependence of the nonlinearity g. 



FIG. 14: (Color online) Asymmetry of the nonlinear RET 
peaks. Shown is the decay rate T(F) of the two peaks (left 
and right) around the respective peak positions F rcs for g — 
(black dashed) and g — —0.3 (solid blue line). 



FIG. 15: Pulsed output from a tilted optical lattice for Vo = 
0.8, F = 1/18 and S — 0,g = 0. Shown is the atomic density 
|^>(x,t)| 2 in a greyscale plot. 



the right one, i.e. the peak decay rate is larger for smaller 
external fields F. This phenomenon becomes even more 
pronounced in the nonlinear case g < 0. The stabiliza- 
tion by an attractive nonlinearity (cf. |13|) is stronger 
for the right peak. 

The peak positions shown in the lower panel of figure 
1 131 vary linearly with g, which can be derived in a pertur- 
bative approach. As discussed above, resonant tunneling 
is observed when a state of a lower ladder get in reso- 
nance with a state in a higher ladder at a different site, 
E a ,n = E a i n i . Here we consider only the states in the 
ground ladder, which are localized in a single potential 
well. First order perturbation theory with respect to the 
linear case g = predicts that their energy is shifted by 
the amount 



0,0 



+ 7V 



g / \4%)\ 4 dx, 



(24) 



where the superscript (0) refers to the linear case g — 0. 
The shift AF of the RET peaks then follows from the 
modified resonance condition 



E 



(0) 



• AE, 



0,71 



E 



(0) 



n'dAF. 



(25) 



Evaluating the integral in equation 124|) and setting n' = 
n — 1 one finds 



0.36 

Z7T 



(26) 



for the peak shift plotted in figure Both peaks are 
shifted equally, so that the distance of the peaks remains 
constant. 

Furthermore the shape of the RET peaks becomes 
asymmetric in the nonlinear case. This is shown in figure 
1141 where we have plotted a magnification of the decay 
rate T(F) around the respective positions of the RET 
peaks for g = —0.3 and g = 0. In comparison to the 
linear case, the peak is bent to the left for an attractive 



nonlinearity. This asymmetry is a general feature of non- 
linear eigenstates in open systems. It is already present 
for the nonlinear two-level system as shown in figure 01 
A similar incline of resonant curves is also important for 
nonlinear resonant transport. The curves can even bend 
over for strong nonlinearities leading to a bistable behav- 
ior as shown in l35l. 



VI. DYNAMICS 

In this section, we discuss the dynamics of an initially 
localized matter wave, e.g., a Gaussian wave packet in the 
tilted double-periodic optical lattice i|19[l . In one of the 
first experiments on the macroscopic dynamics of BECs 
in optical lattices it was shown that such a system shows 
a coherent pulsed output j2^. An explanation in terms 
of truncated Wannier-Stark resonances can be found in 
[36j . The amplitude of the pulsed output is given by the 
decay rate of the Wannier-Stark resonances. 

First of all, we illustrate how the sensitive dependence 
on the phase can be be used to tune a pulsed atom laser. 
As a proof of principle, we just consider the linear (g = 0) 
evolution. In contrast to the previous sections we con- 
sider a weaker potential, Vo = 0.8 so that decay is gener- 
ally stronger. We numerically integrate the Schrodinger 
equation for an initially Gaussian wave packet 



ip(x,t= 0) 



1 



(2 7 r) 1 /4 (7 i/2 



exp(-(x-x ) 2 /4a 2 ) (27) 



with width a = 5ir. FigurelTSlshows the density \ip(x, t)\ 2 
in a greyscale plot for a single periodic lattice (6 — 0). A 
pulsed output forms due to the external field. The pulses 
are accelerated just as classical particles. The first three 
strong pulses emerge from excited ladders. The output 
strength of the other pulses from the ground ladder can 
be controlled to a large extent in the double-periodic case 
by the relative phase <p. As already shown in figure 1101 
the decay rate varies strongly with the relative phase of 
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FIG. 16: Pulsed output from a tilted optical lattice for g — 0. 
Shown is the atomic density for t = 14 Ts for 5 = (top), <5 = 
0.2 and = (middle) and S = 0.2 and = tt/2 (bottom). 



FIG. 17: Fidelity of the output switch in dependence of the 
standard deviation A0 of the phase fluctuations. 



the two lattices. The effect on the pulsed output is shown 
in figure^] where the density is plotted for t = 14Tb for 
three different lattice setups. 

For 6 = 0.2 and = 0, a RET peak of type I is found 
in the second miniladder at F = 1/18. Thus decay from 
this miniladder is strongly enhanced in comparison to 
the single-periodic lattice (5 = 0. Note that the pulsed 
output will stop as soon as the population in the second 
miniladder has decayed. In contrast, the pulsed output is 
strongly suppressed for a relative phase — 7r /2, where 
RET does not play a role for the given field strength. 

In order to measure the performance of this output 
switch more qualitatively, we define the fidelity 

, = -Pout (0o = 0) 

1 Pout(^o=7r/2)' 1 ' 

where P ut measures the integrated density of the pulsed 
output for a certain value 0o of the relative phase of 
the two optical lattices. The output is switched on for 
0o = and it is switched off for 0o = 7r /2. We measure 
the output density P ut at i = 14Tb, where we neglect 
the first three strong pulses as they are due to the initial 
population of excited Bloch bands. Then one has 

/■- 30<i 

P= \i>(x)\ 2 dx. (29) 

J-230d 

For the parameters used in figure 1161 we find a fidelity 
of F = 7.4, i.e., the output for 0o = is enhanced by a 
factor of 7.4 in comparison to 0o = 7r/2. 

However, in a real experimental setup the phase 
will fluctuate around the desired value. This is mainly 
due to mechanical perturbations, thus fluctuations with 
very high frequencies arc unlikely, while fluctuations with 
small frequencies up to some kHz can be controlled by 
an active stabilization. Thus we assume that the power 
spectrum of the phase fluctuations has a maximum at 
intermediate values in the kHz regime. This is compa- 
rable to the Bloch frequency lob — 2tt/Tb, as the Bloch 



period is about one millisecond in a typical experiment 
28]. Exemplarly, we consider fluctuations with a Gaus- 
sian power spectrum with mean ujb and width w.b/4. In 
the following we analyze the pulsed output in dependence 
of the strength of the fluctuations. Figure UTI shows the 
fidelity of the output in dependence of the standard de- 
viation A0 = (((0 — 0o) 2 }) 1 ^ 2 of the fluctuations for the 
same parameters as in figure llfil One observes that the 
fidelity drops to one (no switching effect) for a standard 
deviation of A0 ss 0.67T. A reduction of the fluctuations 
below this value is in principle possible today, however 
only with a great technical effort. As a consequence, an 
output switching seems feasible in double periodic lat- 
tices. 



VII. CONCLUSION AND OUTLOOK 

In the present paper we have studied the interplay be- 
tween decay and a nonlinear mean-field potential describ- 
ing the atom-atom interactions in a dilute Bose-Einstein 
condensate. 

As an illustrative model we have investigated a two- 
level systems with one decaying level, which can be 
treated analytically. In the linear case, one has to dis- 
tinguish two types of level crossings, either the real parts 
anticross while the imaginary parts of the eigenvalues 
cross (type I) or the other way around (type II). Both 
real and imaginary parts are degenerate at the excep- 
tional point, where the bare state decay rate equals the 
coupling strength. A weak nonlinearity does not alter the 
crossing type, however it deforms the levels in a charac- 
teristic manner. For example, it leads to a bending of 
the peaks in the decay rates. Novel nonlinear eigenstates 
emerge for a stronger nonlinearity, where the critical non- 
linearity is decreased in the presence of decay. Looped 
levels appear for type I crossings, while the additional 
eigenstates emerge in a double S-structure for a type 
II crossing. At the exceptional point novel eigenstates 
emerge, even for small nonlinearities. 
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An experimental setup where both decay and nonlin- 
earity play an important role is the dynamics of Bose- 
Einstein condensates in accelerated optical lattices. In 
particular, we have analysed the decay in a double- 
periodic lattice, where a weak period-doubled potential 
is superimposed onto the fundamental lattice. These re- 
sults will be of interest for controlling transport of ultra- 
cold atoms in future and on-going experiments |37| . 

The decay rate in a double-periodic Wannier-Stark sys- 
tem depends sensitively on the system parameters, such 
as the relative amplitudes of the lattices and the relative 
phase, which can be varied over a wide range. In par- 
ticular, the resonant tunneling peaks of the decay rate 
r(F) split up into two subpeaks. Varying the system 
parameters one can tune these peaks and even achieve 
a crossover from a type II crossing to a type I crossing. 
This could be crucial for future experiments since the 
a robust control of the relative phase is hard to realize. 
One can, however, also exploit this sensitive dependence 



in order to implement a fast output switch for a pulsed 
atom laser. A weak nonlinear mean-field potential de- 
scribing the atom-atom interactions in a Bose-Einstein 
condensate of ultracold atoms has two major effects: The 
resonant tunneling peaks are shifted. This shift can lead 
to a stabilization against decay. Furthermore it leads 
to a bending of the peaks as predicted by the nonlinear 
non-hermitian two-level system. 
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